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Abstract. We present new algorithm for growth of non-clustered planar 
graphs by aggregation of cells with given distribution of size and con- 
straint of connectivity k = 3 per node. The emergent graph structures 
are controlled by two parameters — chemical potential of the cell aggre- 
gation and the width of the cell size distribution. We compute several 
statistical properties of these graphs — fractal dimension of the perime- 
ter, distribution of shortest paths between pairs of nodes and topological 
betweenness of nodes and links. We show how these topological prop- 
erties depend on the control parameters of the aggregation process and 
discuss their relevance for the conduction of current in self- assembled 
nanopatterns. 



1 Introduction 

In recent years increased interests in various networks realizations |1I2| revealed 
that several new types of graphs termed structured graphs are more appropriate 
mathematical objects to describe complex network's geometry than traditional 
random graphs 3 . The variety of structures was found to emerge through evo- 
lution processes in which nodes and links are added sequentially according to 
specified rules, in particular, the preferential attachment rules lead to strongly 
inhomogcncous scale-free graphs 1 . In contrast to the evolving networks, which 
comprise a class of causal graphs, the class of homogeneous graphs consists of 
graphs with fixed number of nodes and fluctuating or rewiring links according 
to given rules or certain optimization processes. Complex graph structures may 
emerge in this procedures, especially when certain global or local optimization 
constraints are imposed 0]. 

Planar graphs are special class of graphs that can be embedded in a Euclidean 
plane. A graph is planar iff it does not contain a subdivision of K$ (5-clique) and 
Ka t 3 (minimal non-planar graph with 6 nodes) [3]. Consequently, planar graphs 
fulfill Euler's law: N p + N = E + 1, which is relating the number of nodes N, 
links E and polygons N p . 

In this work we suggest a new method for growing a planar cellular graph 
by attachment of objects — cells (polygons) of length n p , which are chosen from 
a given distribution f(n p ). The polygons are added sequentially in time starting 
from an initial polygon. In addition, we strictly impose the constraint on number 



of links per node k = 3, which is thus fulfilled everywhere in the interior of the 
graph and on some nodes on the graph boundary. The attachment of cells is 
controlled by two parameters — the width of the distribution of cell sizes /i2 and 
the parameter v that plays the role of chemical potential of cell aggregation. In 
the limit of vanishing attachment potential v — ► the growth process resembles 
the one in diffusion- limited aggregation 6 . However, aggregated are spatially 
extended cells of particles rather than single particles. 

Emergent structures of cellular networks are resembling of soap froths [7] or 
patterns of nano-particlcs self- assembled through nonlinear dynamic processes 
8 9 . Typically, a pattern of cells appears when nano-particles are immersed in 
a liquid film, which is then allowed to evaporate until holes of different sizes 
open-up leaving particles in the walls between the holes |9ll(Jj . Generally, the 
structure of the patterns effects the physical processes on them, such as current 
transport . It is therefore important to understand the topology of the ag- 
gregated cellular networks in detail. Here we study the topological properties, 
such as shortest paths between nodes, topological centrality, and fractality of the 
graph's perimeter in different cellular networks obtained by varying the control 
parameters of the aggregation processes. 

2 Cell Aggregation 

The basic idea is to make growing model of planar graph with given distribution 
of cell (polygon) sizes f(n p ) and with degree of nodes 2 < k < 3. 

2.1 Topological constraint 

For this purpose we observe some topological constraints on the distribution 
of cell sizes: (i) f{n p ) is defined for n p > 3, for non-clustered graph we fix 
/(3) = 0; (ii) planar graph obeys Euler's law: N p + N = E + 1. Among these the 
homogeneous plane-filling structures are of special interest |7]. For this class of 
graphs majority of nodes are in the interior of the graph, ie., nodes with degree 
k = 3. Therefor 3N » 2E and Euler's law becomes 

6N p = 2E + 6. (1) 

For large system with distribution of cell sizes f(n p ) we have 

N = Np J2^^, E = N p J2 Vl ^- (2) 

Tip Tip 

Substituting second relation into QJ: 6N P — N p J2 n n pf( n p) + 6, then for large 
N p 3> 1 one can find that the average cell size is equal to six 



( n p) = ^2 n pf( n p) = 6. 
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We use lognormal distribution of polygon size 
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which is most often found in experiments |2| ■ Using the condition the number 
of independent parameters in 10} is reduced 



where second central moment \ii remains as the control parameter in our case. 
2.2 Model 

Starting from an initial cell, at each time step a new cell with size taken from 
the distribution f(n p ) is attached to the graph boundary of the graph according 
to the following rules: 

1 . In the aggregation only nodes with degree 3 can be closed inside the graph; 

2. New link can be assigned only to nodes with degree 2; 

3. Potential nesting place is searched as an array of nodes on graph boundary 
with degree 3 limited with two nodes of degree 2. The nested part of the cell 
is identified with the nodes of the nesting string. Therefore, number of extra 
nodes to be added is n = n p — I, where I is length of the nesting string. We 
select the nesting place with probability p ~ e^ 1 " 1 , where the parameter v 
plays the role of the chemical potential for addition of new nodes. 
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Fig. 1. Two possible types of cell-aggregated planar graphs: open and closed 
structures, obtained by lognormal distribution with /12 = 1.0, and aggregation 
potential v — 5.0. 



We implemented this algorithm in C+ + according to following steps: 

Initial graph: one polygon of size np taken from f (np) 
For i=2 to Np 

np = next random from distribution f(np) 

If (there is no nodes on graph boundary with degree 2) exit(l) 
For all j=(periphery node with degree 2) 

d=distance to the next node on graph boundary with degree 2 

Number of new nodes n = np-d-1 

lf(n>0) p(j)=exp(-nu*n) 
End of loop j 

If (there is no j with n>0) exit (2) 
Normalize p(j) 

j = next random from distribution p(j) 

Add new polygon with size np linked with 

node j and next node on graph boundary with degree 2 
End of loop i 
exit(O) 

Depending on model parameters of the growth process and it stochasticity 
three possible exit cases are: 

— exit(O) - Open structure (planar graph with N p polygons); 

— exit(l) - Closed structure (after some number of step there are no more nodes 
of degree 2 and structure stops to grow, no nesting places of any size); 

— exit(2) - No nesting place available for current cell. In this case one can take 
next cell, which in turn perturbs the actual distribution. 

We never experienced the exit (2) situation for the range of parameters [ii G 
[0.5, 2.0], v e [0, 5] and N p — 1000 in huge number of samples. Two examples of 
the emergent open and closed structures are shown on Fig. ^ More examples of 
cellular networks are shown in Fig. [3 for varying parameters /i2 and v.. 

3 Fractal Dimension of Network Perimeter 

During the nesting growing process in one step number of nodes N increases 
by n < n p . Number of added nodes n at each step depends on the cell size and 
length of the nesting string. Therefore, N w nN p where k — (n) is average growth 
rate. For open structures (see Fig. boundary of the graph becomes fractal, 
depending on the control parameters. In fact, length of the graph boundary L 
grows as a power of the number of cells N p (or network size N) with fractal 
dimension D defined by L ~ N p . 

In Fig.|3we show how number of nodes on the graph boundary increases with 
N p . Each point is averaged over 10 emergent growing networks. The dimension 
D is in the range | < D < 1, when D = 1 correspond to structures of high 
fractality, that is obtained for small values of the parameter v. D = | correspond 
to planar " circle like" structures with reduced fractality. For > we observe a 
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Fig. 2. Cell- aggregated planar graphs with lognormal distribution of polygon 
size for various values of width /12 and cell aggregation potential v. Width of 
lines represent topological betweenness (centrality) of links calculated in Section 
[PI 
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Fig. 3. Scaling of the network perimeter for (a) ^2 = hexagons only (b) /K2 
2.0 



continuous crossover between these two limits (see Fig.[3t> for // 2 = 2). However, 
in structures with homogeneous cell distribution (p,2 = - hexagons only) a 
sharp transition seems to occur at v c w 1.5. 



4 Shortest Paths and Centrality on Cellular Networks 

In this section we consider global topological properties of the cell-aggregated 
planar graphs and their dependence on the control parameters //2 and v. 

4.1 Shortest Paths 

Shortest path between two nodes is defined as path along the smallest number of 
intermediate links 3 . We implemented an algorithm for counting shortest paths 
of Dijksta type I n Fig. 0] we show distribution of lengths of shortest paths 
between all pairs of nodes on network. All networks are for fixed n% = 1.0 and 
approximately of the same size N sa 1000 nodes. Each point in Fig.^Jis averaged 
over 100 sample networks. We found similar results for other fi2 values. 
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Fig. 4. Distribution of lengths of shortest paths on networks for fixed \ii = 1.0 
and various values of v. Insert: The distribution in the case v = is shown on 
linear scale. Solid line: Gaussian with Iq = 23 and a = 8.37. 



All these networks have similar topology at local level, because the number of 
links at all interior nodes is constant k = 3. Therefore, distributions of shortest 
distances at small scale are similar for all values of parameter v. Differences in 
global topology appear on large scale for lengths larger then peak value lo ~ 25, 
which manifest in occurrence of additional peaks (see Fig. The probability 
of long paths increases for larger values of the parameter v. Whereas in the 
limit v = the distribution of length of shortest paths on large scale can be 
approximated with a normal distribution (inset on Fig.^J. 



4.2 Centrality Measures 



Bctweenness centrality of a node in network is defined by 

Cb(v) = £ ^ (6) 

where cr s t is total number of shortest paths between nodes s and t, and cr s t(v) 
is number of these paths that node v lies on. Betweenness of links is defined in 
analogous way. In our algorithm, we record number of shortest paths through 
each node and through each link on a network. In Fig. [5] we show distributions 
of betweenness of nodes and links averaged over 100 sample networks with fixed 
fi2 = 1-0 and size N w 1000 nodes. For inhomogeneous networks, distributions 




Fig. 5. Distribution of betweenness (a) of nodes (b) of links. 



of these two betweenness measures can be substantially different. In our case, 
however, they are similar because all interior nodes have fixed degree k = 3. 
We find that distributions at small scale strongly depend on parameter v, which 
results from the graph boundary. Similar conclusions hold for other fj,2 values. 

In Fig. [21 thick lines represent links with large betweenness. For this type of 
networks the strongest lines, which make the skeleton of the graph, are connect- 
ing the nodes with largest centrality measure. 



5 Conclusions 



We have introduced a new algorithm for growth of graphs by aggregation process 
of extended objects - polygons with size distribution. Depending on aggregation 
conditions, which are determined by two parameters controlling the distribu- 
tion width and attachment potential, we can get a wide spectrum of emergent 
structures. In this paper we presented some results for the case of lognormal 
distribution of cells and additional constraints, leading to the emergent non- 
clustered planar graphs with a constant node connectivity. The algorithm works 
for variety of cell distributions and constraints, that may result in diverse opened 
(fractal) or closed structures. For instance, for a special set of parameters we can 
get closed structures of Cqq type. 

Wc measured several topological properties of these networks in quantitative 
details — fractality of the graph boundary, shortest paths, and betweenness cen- 
trality. This properties are important for some dynamical process on networks 
such as electrical conductivity JT] via single electron tunneling |9I14| . 
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